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ABSTRACT 


During  FY  79,  the  NBS  Thermophysical  Properties  Division  program  of  research  for 
superconducting  power  transmission  line  (SPTL)  development  focused  on  three  tasks: 

1)  Analytical,  numerical,  and  experimental  modeling  of  counter-flow  SPTL  cool-down. 

2)  Measurement  of  radial  permeability  of  a SPTL  cable  to  gas  flow. 

3)  Experimental  evaluation  of  thermal  flux  meters  as  a possible  technique  for  deter- 

mining enclosure  heat  leak. 

We  have  completed  the  first  phase  of  our  investigation  of  SPTL  cool-down  in  which  the 
aim  was  to  develop  a fundamental  understanding  of  the  counter-flow  method  of  cool-down. 

The  excellent  agreement  of  our  analytical,  numerical,  and  experimental  results  gives  us 
confidence  that  we  can  accurately  model  the  cool-down  of  full  scale  SPTL's. 

The  cable  permeability  measurements  indicate  that  any  significant  rupture  of  the  lead 
gas  barrier  of  the  cable  will  lead  to  an  unacceptably  high  leak  rate  from  the  pressurized 
core. 

Evaluation  of  thermal  flux  meters  continues,  but  preliminary  results  indicate  large 
effects  due  to  seasonal  variations  in  the  soil  heat  flux  and  percolation  water  through  the 
soil. 

During  FY  80  we  shall  concentrate  on  developing  cool-down  strategies  for  full  scale  SPTL's 
using  the  detailed  computer  code  previously  developed  at  NBS.  Field  evaluation  of  thermal 
flux  meters  will  also  continue. 

KEY  words:  Cable  cool-down,  cable  permeability,  counter-flow  cool-down,  heat  flux  meters, 

liquid  helium,  permeability  of  composites,  superconducting  power  transmission. 
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1.0  TRANSMISSION  LINE  COOL-DOWN 
(D.  E.  Daney,  P.  R.  Ludtke  and  M.  C.  Jones) 


1.1  Introduction 

The  time  required  to  cool  down  a superconducting  (or  cryoresi sti ve)  power  transmission 
line  is  critical  to  the  success  of  these  lines  since  excessive  cool-down  and  warm-up  times 
would  result  in  unacceptably  long  interruptions  of  service  when  repairs  are  required.  For 
example,  Jones  [1]  has  calculated  a cool-down  time  of  20  days  for  a 15  km  line  cooled  by  a 
two  stream  counterflow  method  (thermally-coupled  return).  Since  the  warm-up  process  is 
simply  the  inverse  of  the  cool-down  process  (at  least  for  the  constant  property  case),  a 
service  interruption  of  six  weeks  or  more  could  be  expected  each  time  internal  repair  or 
maintenance  is  required  for  a SPTL  of  this  configuration. 

If  we  understand  the  exact  nature  of  the  cool-down  process,  then  we  may  alter  the 
design  or  operating  conditions  to  give  acceptable  cool-down  times.  For  example,  one 
solution  to  excessive  cool-down  times  is  to  reduce  the  length  between  refrigeration  sta- 
tions since,  as  we  demonstrate  in  this  report,  the  cool-down  time  is  proportional  to  the 
square  (constant  mass  flow  rate)  or  cube  (constant  pressure  drop)  of  the  distance  between 
refrigeration  stations  for  counter  flow  cooling.  As  is  often  the  case  in  engineering 
designs,  the  design  of  the  SPTL  cooling  configuration  may  be  governed  by  transient  rather 
than  steady-state  considerations. 

Two  arrangements  have  been  proposed  for  cooling  cryogenic  power  transmission  lines. 

The  more  usual  arrangement.  Fig.  lb,  has  the  return  stream  thermally  isolated  from  the  go 
stream  (thermally-uncoupled  return).  Because  this  method  requires  an  expensive  vacuum- 
insulated  return  line.  Dean  and  Jensen  [2]  proposed  a system  with  thermal  contact  between 
the  go  and  return  streams.  Fig.  la,  (thermally-coupled  return)  for  the  Brookhaven  National 
Laboratory  (BNL)  line.  This  same  general  scheme  has  also  been  adopted  in  the  Los  Alamos 
Scientific  Laboratory  (LASL)  proposal. 

In  this  latter  arrangement  supercritical  helium  flows  down  the  hollow  core  of  the 
cable,  is  expanded  to  a lower  pressure  (and  temperature)  at  the  far  end,  and  then  is 
returned  in  the  space  outside  the  cable  within  the  enclosure.  The  wrapped  dielectric 
provides  weak  local  thermal  contact  between  go  and  return  streams.  Integrated  over  a 
length  of  10  km  or  15  km,  however,  this  configuration  produces  a rather  good  counterflow 
heat  exchanger,  N^^  ~ 40. 

The  cool-down  characteristics  of  these  two  arrangements  are  markedly  different  as 
indicated  in  Figure  1.  In  the  case  of  the  thermally-uncoupled  return,  a cool-down  wave 
propagates  along  the  length  of  the  line,  and  cool-down  is  complete  when  this  wave  has 
passed  the  far  end.  In  contrast,  in  the  thermally-coupled  case  (counterflow)  a cool-down 
(temperature)  wave  developes  along  the  whole  length  of  the  line,  and  cool-down  is  not 
complete  until  the  amplitude  of  this  wave  decays  to  zero.  The  rate  of  decay  is  propor- 
tional to  the  temperature  difference  between  the  go  and  return  streams  (i.e.,  to  the  rate 
of  refrigeration)  which  in  turn  is  proportional  to  the  thermal  resistance  between  the  go 
and  return  streams. 


INSULATED  RETURN 


a) 


Thermally-coupled  return  both  ends  cooled  b)  Thermally-uncoupled  returns 

Figure  1.1  Configuration  and  cool-down  characteristics  of  cooling  schemes. 


Prior  to  the  proposal  of  Dean  and  Jensen  [2]  for  the  BNL  line,  Wilkinson  [3],  Garwin 
and  Matisoo  [4],  and  Norris  and  Swift  [5]  discussed  the  thermally-coupled  return  arrange- 
ment. Flemming  [6],  Edney,  et  al.  [7]  and  Morgan  and  Jensen  [8]  analyzed  the  steady  state 
operating  behavior  of  this  arrangement.  For  constant  properties  and  symmetry  in  the  heat 
transfer  and  boundary  conditions,  a parabolic  temperature  distribution  results.  When  the 
thermal  contact  is  good  between  the  go  and  return  helium  streams,  a large  temperature 
maximum  will  occur  at  the  midpoint  of  the  line.  This  occurs  because  as  the  thermal  contact 
between  streams  increases,  the  net  enthalpy  flux  of  coolant  across  a section  of  the  line 
decreases.  The  same  phenomena  occurs  during  cool-down,  except  in  this  case,  the  heat  load 
consists  primarily  of  the  line  heat  capacity. 

The  object  of  this  study  has  been  to  investigate  the  cool-down  behavior  of  supercon- 
ducting power  transmission  lines  which  use  the  thermally-coupled  return  arrangement.  Our 
original  commission  was  for  detailed  numerical  modeling  of  cool-down  of  the  BNL  and  LASL 
lines  (reported  previously  by  Jones  [1])  followed  by  experimental  modeling  of  this  thermally- 
coupled  return  arrangement.  We  have  expanded  the  scope  of  our  investigation  somewhat  so 
the  second  phase  of  the  work,  which  is  reported  here,  includes: 

1)  An  approximate  analytical  solution  in  terms  of  dimensionless  variables  for  the 
constant  property,  constant  mass  flow  rate  case. 

2)  Numerical  solution  in  terms  of  dimensionless  variables  for  the  constant  property, 
constant  mass  flow  rate  case. 

3)  Experimental  study  of  the  thermally-coupled  return  arrangement  for: 

a)  cool-down  to  76  K 

b)  cool-down  to  6 K 

c)  warm-up  from  76  K and  6 K 

4)  Detailed  numerical  modeling  of  the  experiments  for  experimental  confirmation  of 
the  computer  program  developed  by  Jones  for  modeling  of  the  BNL  and  LASL  lines. 

The  experimental,  analytical,  and  numerical  results  agree  well,  so  that  we  may  now 
predict  transmission  line  cool-down  times  with  good  accuracy. 
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1.2  Analysis 

For  analysis  of  the  thermally-coupled  return  (counterflow)  arrangement,  we  choose  the 
simple,  one-dimensional  model  shown  in  figure  1.2. 


Figure  1.2  Simple  model  of  transmission  line. 

The  solid  is  lumped  into  a single  wall  separating  the  go  and  return  stream,  and 
constant  properties  and  mass  flow  rate  are  assumed.  The  thermal  resistance  is  concen- 
trated at  the  wall-fluid  interface. 

The  time  dependent  energy  equations  for  the  problem  (obtained  from  the  first  law  and 
continuity  equation  in  the  standard  manner  and  neglecting  axial  conduction  and  any  time 
dependence  of  the  pressure)  are: 

for  the  go  stream 

9T  9T 

^ (C  pA)  dx  = - mC  ^ dx  - (hP)  dx  (T  - T ) 1.1 

9t  ^ p^  '^u  p 9x  ^ '^u  ^ u w^ 

for  the  return  stream 

9T  9T 

(C  pA)  dx  = iTiC  dx  - (hP)  dx  (T  - T ) 1.2 

9t  ^ p^  p 9x  ' V w^ 

and  for  the  wal 1 
9T 

^ (pAC)  dx  = (hP)  dx  (T  - T ) - (hP)  dx(T  - T ) 1.3  ' 

9t  '^w  ^ -^u  u w'^  ^ ^ w u^ 

If  we  set  (hP)^  = (hP)jj  and  A^  = A^,  then  equations  1.1  through  1.3,  expressed  in  terms 
of  dimensionless  variables  (defined  in  Section  1.7),  become 
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u 


1.4 


1 3u  _ 9u 
b 9t  9y 


1 ^ 
b 9x 


9y 


— + w - V 


1.5 


9w 


U + V 


2w 


1.6 


with  initial  conditions: 


u = V = w = 1 for  all  y at  i = 0 


1.7 


and  boundary  conditions: 


u = 0 at  y = 0 
u = V at  y = L 


for  T > 0 


1.8 


Thus,  we  consider  the  response  of  a transmission  line,  initially  at  uniform  temperature, 
to  a step  change  in  the  coolant  inlet  temperature  at  the  near  end. 

The  second  boundary  condition  is  for  cooling  from  one  end  only  (single-ended  cool- 
down). In  this  case  the  coolant  streams  are  short  circuited  at  the  far  end  (x  = £)  so 
that  the  go  stream  is  connected  directly  to  the  return  stream,  resulting  in  T^  = T^  at 
X = Jl.  If  both  ends  are  cooled  (i.e.,  T^  = v = 0 at  y = L)  the  temperature  distribution  is 
symmetrical  about  the  midpoint,  and  the  solution  for  a half  length  of  this  double-ended 
cool-down  corresponds  to  the  solution  for  single-ended  cool-down.  Thus  we  need  only 
obtain  the  solution  for  single-ended  cool-down  and  the  following  analysis  applies  to  this 
case.  At  the  conclusion  of  the  analysis,  transformations  for  double-ended  cooling  are 
given. 

Numerical  integration  of  equations  (1.4)  through  (1.8)  was  accomplished  using  the 
computer  program  PDECOL  developed  by  Madsen  and  Sincovec  [9]. 

We  are  unaware  of  an  analytical  solution  to  this  problem.  We  obtain  an  approximate 
solution  by  neglecting  the  time  dependent  terms  in  equations  (1.4)  and  (1.5)  as  has  been 
done  in  analysis  of  the  regenerator  problem  [10].  That  is,  we  assume  that  b (the  heat 
capacity  of  the  solid  relative  to  the  heat  capacity  of  the  gas)  is  large.  At  ambient 
temperature  this  is  a good  approximation;  at  helium  temperature  it  is  not.  However, 
we  are  able  to  reintroduce  the  dependence  on  b to  give  cool-down  solutions  good  for  a 
wide  range  of  heat  capacity  ratios. 

Before  giving  the  details  of  the  solution  we  will  first  outline  the  method.  We 
begin  by  assuming  a solution  of  the  form 


u = F(x)  • G(y)  . 1.9 

Based  on  our  experimental  and  numerical  results,  we  further  assume  a sinusoidal  form  for 
G(y)  so  that  equations  1.4  thru  1.8  may  then  be  solved  to  give 
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V = F(t)  • H(y) 


1.10 


w = F(t)  • I(y)  . 1. 11 

where  H(y)  and  I(y)  are  known  and  F(t)  remains  to  be  determined.  With  the  spatial  temperature 
di strubutions  known  as  functions  of  the  assumed  spatial  distribution  G(y),  we  may  now 
integrate  the  first  law  of  thermodynamics  for  steady  flow 


mC  (v 
P 0 


u )dt  = -(pAC) 

0 ^ w 


(1-f) 


a 

[w(x,t+dt)-w(x,t)] 


dx 


1.12 


to  obtain  the  time-dependent  part  of  the  solution  F(t).  The  term  2/b  accounts  for  the 
heat  capacity  of  the  two  gas  streams.  The  key  to  the  method  is  to  obtain  an  expression  for 
- u^  (which  gives  the  rate  of  refrigeration)  for  use  in  eguation  1.12. 

From  our  experimental  and  numerical  work,  we  observe  that  counterflow  cool-down  pro- 
gresses in  two  stages.  During  the  first  and  relatively  brief  stage,  a cool-down  wave  propagates 
towards  the  far  end  (y  = L)  with  u^  undiminished,  figure  1.3. 


y 


Figure  1.3  Assumed  temperature  profile  during  first  stage  of  cool-down. 


During  the  second  and  longer  stage,  the  cool-down  wave  extends  the  full  length  of  the  line 
and  the  amplitude  of  the  wave  undergoes  an  exponential  decay,  figure  1.4.  The  instan- 
taneous rate  of  decay  depends  on  the  rate  of  refrigeration  which  is  proportional  to  the 
near  end  fluid  temperature  difference,  v^  - u^. 
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y 

Figure  1.4  Assumed  temperature  profile  during  second  stage  of  cool-down. 


Returning  now  to  the  details  of  the  solution,  we  approximate  the  go-stream  temperature 
distribution  by: 

u = u^  sin  I ^ 1. 13 

where  u = u (t).  Eg  (1.4)  gives: 

w = u^  (sin  ^ ^ cos  I 1.14 

Intergration  of  equation  1.5  gives: 


2 u. 


V = 


1 + 


2 T 2 

/l\l^  1 ^ y L ^ 

(-^)  exp  (y-L)  + 2 1 " (-Jl^  ? i T 


2 L 


^ y 

2L  2 L 


1.15 


Evaluating  v at  y = 0,  we  obtain: 


V 

0 


7T 

L 


1 + 2L  exp  (-L) 
''  2L^ 


1.16 


6 


which,  noting  that  e 0,  is  the  normalized  temperature  difference  (v^  - u^)  between  the 
go  and  return  streams  at  the  cooled  end  of  the  transmission  line.  For  large  L,  equation 
(1.16)  reduces  to: 


V 

0 


n 

L 


1. 17 


We  will  now  use  this  result  in  equation  1.12  to  obtain  the  time-dependent  part  of  the 
solution.  During  the  first  stage  of  cool-down,  we  assume  the  cool-down  wave  progresses 
towards  the  end  of  the  line,  as  shown  ’in  figure  1.3,  such  that: 

w = sin  ^ ^ for  y<k  1.18 


and 


w = 1 for  A<y<L 


1.19 


where  A = A(t)  is  the  position  of  the  cool-down  wave  front. 


Although  (1.18)  does  not  strictly  agree  with  (1.13),  its  use  simplifies  evaluation  of 
the  integral  in  (1.12).  The  error  in  the  total  cool-down  time  introduced  by  this  approxi- 
mation is  small  because  the  duration  of  the  first  stage  is  short  compared  to  that  of  the 
second  stage. 

Substituting  (1.17)  and  (1.18)  into  equation  (1.12)  and  noting  that  1 during  the 

first  stage  gives 


/ 


sin  5 ^^dy  + (L  - Ag) 


/ 


sin  ^ ^^dy  - (L  - AJ 


1.20 


or 


S dl  = (1  - f)  (1  + f)  dA 


1.21 


and  finally. 


^1  = 


1.22 


as  the  dimensionless  time  for  the  first  phase  of  the  cool-down. 

The  second  stage  cool-down  time  is  obtained  in  an  analogous  fashion.  The  first  law 
(including  the  heat  leak  Q and  the  heat  capacity  of  the  gas)  gives: 


(- s 


+ 


dt  = - (pAC)^  (1  + |) 


/ 


Jl 

[w(x,t+dt)  - w(x,t)] 


dx 


1.23 
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Replacing  w and  by  equations  (1.14)  and  (1.17)  and  transforming  to  dimensionless  vari- 
ables gives: 


ITlC 


(T.  - T^) 


^ n 
L 


dt2  = 


du„  (1 


f) 


/ • ^ y 

(sin  2 f 


n 

2^  cos 


5 y.) 

2 L-' 


dy 


1.24 


Evaluating  the  integral  and  rearrangement  yields: 


CSI 

1 

II 

Q 

7T 

di2 

2L^  (1  + II) 

" b 

mCpUT.  - Tp^ 

_ 2(1*  ir) 

1.25 


which  has  the  solution: 


1 - 


_QL 


nftiCpdi  • 


exp 


71^  "^2 


2L^(1  + ^ b 


1.26 


niiiC  (T.  - T,)„ 
p^  1 f^SL 


for  the  second  stage  of  cool-down. 

To  summarize  the  cool-down  behavior  of  the  thermally-coupled  arrangement,  we  find 
that  following  an  initial  lag  given  by  (1.22),  the  normalized  end  temperature  u^  (note 
that  ~ “ v^)  undergoes  an  exponential  decay  from  u^  = 1 to  u^  = QL/;n[iiiCp(T^.  - T^)^. 

The  term  QL/TtiiiCp  is  the  rise  in  the  final  steady  state  temperature  of  the  end  of  the  line 
due  to  heat  input.  It  is  the  maximum  temperature  in  the  line  and,  it  has  a value  L/n 
times  the  temperature  rise  that  would  occur  with  no  thermal  coupling  between  the  two 
streams  (i.e.,  the  temperature  rise  for  the  go  stream  only  - single  ended  cooling). 

We  may  compare  this  estimate  of  the  maximum  steady  state  temperature  with  the  exact 
solution  for  the  steady  state  constant  property  case  with  heat  leak  [7],  which  is 


(1. 


L)^ 


(-Q-) 


4L 


1.27 


for  L > 1 with  equal  heat  inputs  to  the  go  and  return  streams.  Evaluating  the  expres- 
sions, we  obtain: 


L 

(1  + L)^ 

L 

4L 

71 

1 

1.00 

0.32 

n 

1.36 

1.00 

10 

3.03 

3.18 

20 

5.  51 

6.37 

50 

13.01 

15.90 

00 

L/4 

L/n 

The  agreement  of  the  sinusoidal  approximation  with  the  exact  solution  is  adequate  for  our 
purposes,  because  in  most  cases  of  practical  interest,  the  temperature  rise  due  to  heat 
leak  will  be  small.  Our  numerical  results  indicate  that  the  sinusoidal  approximation 
should  be  limited  to  cases  where  L > 5. 


The  cool-down  constant  time  given  by  (1.26)  is: 


2_ 

nZ 


5_  ) (2  + -)L^ 
2L  ’ 


1.28 


for  the  second  stage  of  cool  down.  The  constant  N is  2.3  for  90  percent  cool-down  and  4.6 

for  99  percent  cool-down  (i.e.,  90  or  99  percent  of  u - QL/niiiC  ). 

a p 

The  preceeding  results  apply  to  cooling  from  one  end  (single-ended  cooling).  They 
may  also  be  applied  to  cooling  from  both  ends  (double-ended  cooling)  by  noting  that  the 
appropriate  length  in  the  preceeding  is  ij  the  total  length  because  of  the  symmetry  of  the 
solution  about  the  mid-point. 

When  considering  an  actual  line  or  experimental  measurements,  it  is  more  convenient 
to  use  the  overall  heat  transfer  coefficient  U.  The  corresponding  dimensionless  length 
is  then  the  number  of  heat  exchanger  transfer  units  where  N^^=UA/mCp.  Noting  that 

U = h/2,  the  following  transformations  apply: 

one  end  cooled  both  ends  cooled 


t' 

= 

1/2  T 

1/2  X 

^.u 

= 

1/2  L 

L 

t'/N^ 

tu 

2 t/L^ 

1/2  t/L^ 

The  total  cool-down  time  obtained  adding  (1.22)  and  (1.28)  becomes 


t'  = 1.048 
cd 


tu 


(1-^)  (1-f) 

tu 


1.29 


for  single-ended  cool-down,  and: 

t'  = 0.262  N?  (1  + J^)  (1  + h 
cd  tu  N.  ^ ' b 

tu 


1.30 


for  doubled-ended  cool-down  (90  percent-coo 

We  note  from  1.29  and  1.30  that  for  la 
2 

t'  ,/N.  is  nearly  constant.  If  we  choose  t 
cd  tu 

of  w vs.  t'/N.  should  coalesce  the  results 
tu 

dependence  on 

The  first  order  dependence  of  the  coo' 
directly  from  1.29  or  1.30.  For  large 

t 

cd  m 

For  cases  where  the  total  pressure  drop  acros 


P Ap 


for  an  ideal  gas  coolant.  Expression  1.32  ma; 
2 

that  m ~ p Ap/£  for  an  ideal  gas.  Because  1 
rigorous  proof  of  1.32  requires  integration  o 
the  appendix  (section  1.9). 


down  for  both  cases  i.e.,  w=0.1). 
ge  (say  greater  than  10)  the  parameter 
as  a correlating  parameter,  then  a plot 
to  a narrow  band  of  curves  having  only  a weak 

down  time  on  transmission  line  length  follows 

1.31 

the  transmission  line  is  fixed,  we  obtain 

1.  32 

be  obtained  directly  from  1.31  by  noting 
31  was  derived  for  constant  lii,  however,  a 
1.12  with  iTi  = m(t).  This  proof  is  given  in 
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1.3  Experimental 

The  analysis  shows  that  an  experimental  model  of  a counterflow  (thermally-coupled 
return)  cooled  SPTL  should  be  scaled  so  as  to  have  the  same  values  of  and  b (ratio  of 
cable  to  coolant  heat  capacity)  as  the  full  scale  line.  Typical  ranges  of  these  parameters 

3 

for  a SPTL  are  10  to  100  with  b on  the  order  of  10  (ambient  temperature  values). 

The  construction  and  instrumentation  of  a concentric-channel-test-section  present  a 
number  of  difficulties.  In  a SPTL,  the  go  and  return  streams  are  separated  by  the  dielec- 
tric annulus  (figure  1.1)  which  is  several  centimeters  thick.  Instead,  we  chose  the  con- 
struction shown  in  figure  1.5  wherein  the  go  and  return  streams  are  separated  by  a web 
which  gives  a controlled  thermal  resistance  between  the  two  streams.  Although  geometrically 
different  than  a SPTL,  the  characteristic  dimensionless  cool  down  parameters  are  similar:  b 
= 450  (at  ambient  temperature)  and  ranges  from  10  to  22. 

The  16  m test  section  (figures  1.5  and  1.6)  was  swaged  to  the  two  passage  shape  from  a 
19  mm  X .5  mm  wall,  type  304  stainless  steel  cylindrical  tubing.  The  two  passages,  6.4  mm 
and  4.8  mm  ID,  are  separated  by  a 9 mm  x 1 mm  web  which  has  a continuous  seam  weld  for 
hermetic  separation  of  the  channels.  A 0.05  mm  electro  deposited  layer  of  copper  enhances 
the  thermal  conductance  between  the  passages.  For  compactness,  the  test  section  is  wound 
in  a 37  cm  diameter  coil  (Figure  1.6).  Thermal  insulation  is  provided  by  a high  vacuum’ 
environment,  gold  plating,  10  layers  of  aluminized  film,  and  a radiation  shield  which  may 
be  controlled  down  to  90  K. 

Helium  gas  is  supplied  to  the  test  section  via  a counterflow  heat  exchanger,  a LN2 
bath,  a second-counter  flow  heat  exchanger,  and  a LHe  bath.  The  valve  arrangement  permits 
either  single  stream  cool-down  or  warm-ups,  as  well  as  counter  flow  cool-down.  For  the 
experiments  in  which  the  final  temperature  was  77  K,  liquid  nitrogen  is  used  in  place  of 
the  liquid  helium  shown  in  figure  1.5.  Single-ended  cool-down  experiments  are  run  by 
connecting  the  go  stream  discharge  directly  to  the  return  stream  inlet. 

Fast  thermal  response  of  the  go  stream  and  return  stream  inlets  is  assured  by  close 

coupling  the  go  stream  control  valve  and  return  stream  liquid  reservoir  to  the  test  section 

with  thin-walled  stainless  steel  tubing.  Typical  response  times  of  the  inlet  conditions 

are  2 to  4 percent  of  the  test  section  cool-down  times.  Figure  1.7  compares  the  inlet 

temperature  decay  rate  chosen  for  the  numerical  calculations  to  the  experimental  inlet 

temperature  decay  rate  for  three  different  experiments.  As  discussed  at  the  end  of 

2 

Section  1.2  we  use  the  dimensionless  parameter  (dimensionless  time  divided  by 

dimensionless  length  squared)  because  use  of  this  parameter  provides  a common  base  of 
comparison  for  widely  varying  conditions.  The  agreement  of  the  experimental  inlet  temper- 
ature decay  rates  with  the  decay  rate  assumed  in  the  numerical  calculations  is  adequate. 

The  exact  value  is  unimportant  so  long  as  the  inlet  temperature  decay  time  is  short 
compared  to  the  cool-down  time. 

1.3.1  Instrumentation 

Temperatures  are  measured  with  0.13  mm  diameter  gold  (0.07  at.  percent  iron)  vs. 
chromel  thermocouples  referenced  to  a liquid  nitrogen  bath.  Those  for  the  test  section 
are  mounted  on  the  web  equidistant  from  the  two  flow  channels.  Differential  thermocouples 
are  used  to  measure  the  go  and  return  stream  temperature  differences  at  each  end.  Pres- 
sures are  measured  with  variable  reluctance  diaphram  transducers  operating  at  ambient 
temperature.  The  return  stream  exit  mass  flow  rate  is  measured  with  an  ambient  tempera- 
ture thermal -type  flow  meter.  The  total  uncertainty  in  the  measurements  is  estimated 
to  be:  temperature,  0.5  K;  differential  temperature,  0.05  K;  pressure  0.5  percent; 
and  mass  flow  rate,  0.5  percent. 
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Figure  1.5  Test  section  and  flow  schematic  for  cool-down  experiments. 


Figure  1.6  Dual  passage  test  section,  instrumented  and  wraped  with  ten  layers 
aluminized  film. 


Figure  1.7  Decay  of  fluid  inlet  temperature. 


Figure  1.8  Single  stream  cool-down  experiment  for  evaluation  of  the  test  section  heat 
capacity.  The  curves  are  to  guide  the  eye  through  the  experimental  data. 
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Data  are  collected  with  a mini  computer  which  gives  "real  time"  conversion  of  all 
instrument  channels  to  physical  parameters  (P,T,!t). 


1.3.2  Procedure 

The  study  requires  three  types  of  experiments.  In  addition  to  the  counterflow  cool- 
down and  warm-up  runs,  we  also  measure  the  test  section  heat  capacity  and 

The  test  section  heat  capacity  is  determined  from  single  stream  cool-down  experiments 
using  the  expression: 


(^)  S,  = 


■ Tout’  *TpClt 


1.33 


The  experiment  begins  with  the  test  section  at  an  initial  uniform  temperature  T^. . Single 
stream  cool-down  is  initiated  by  opening  V2  (figure  1.5)  and  closing  VI  with  V4  closed  and 
V5  open.  Figure  1.8  shows  the  results  of  such  an  experiment.  The  integral  in  (1.33)  is 
evaluated  up  to  the  time  when  falls  to  4 or  5 percent  of  its  maximum  value,  i.e., 

until  the  cool-down  wave  has  essentially  passed  through  the  test  section.  The  contribution 
of  heat  leak  to  the  error  in  heat  capacity  is  estimated  to  be  less  than  five  percent. 

The  thermal  coupling  between  the  go  and  return  streams  is  determined  from  the  steady 
state  operation  of  the  test  section  as  a counter-flow  heat  exchanger,  i.e.,  the  go  stream 
inlet  is  held  at  a constant  temperature  near  ambient,  and  the  return  stream  inlet  is  held 
at  76  K or  4 K.  The  and  UP£  are  then  determined  from  the  heat  exchanger  effectiveness: 


N 


^ UP£  ^ £ 

tu  m^  T-e 
P 


1.34 


Figure  1.9  shows  the  measured  test  section  as  a function  of  1/rti.  The  time  required  to 
achieve  steady-state  is  approximately  equal  to  the  cool  down  time  for  the  cool  down  experi- 
ments. 

The  heat  leak  to  the  test  section  during  the  cool-down  experiments  is  evaluated  by 
solving  equation  (1.23)  for  Q.  Thus  the  heat  leak  is  the  difference  between  the  rate  of 
refrigeration  and  the  rate  of  cooling.  Figure  1.10  gives  the  heat  leak  evaluated  by  this 
method  for  three  typical  experiments.  The  complex  time  dependence  of  the  heat  leak  results 
from  the  experimental  arrangement.  Although  the  test  section  and  radiation  shield  are 
cooled  down  together,  the  temperature  of  the  test  section  varies  in  both  position  and 
time.  Consequently,  the  radiation  shield  temperature  matches  the  test  section  temperature 
at  only  a single  position  or  pair  of  positions,  and  the  heat  leak  varies  with  both  time 
and  position. 


Figure  1.9  Test  section  vs.  1/iii  as  determined  from  steady  state  counter  flow 
experiments.  The  curve  is  to  guide  the  eye  through  the  experimental  data. 


TIME,  s 

Figure  1.10  Measured  heat  leak  during  three  typical  cool-down  experiments. 
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1.4  Results  and  Discussion 


Table  1 summarizes  the  experiments.  Abortion  of  two  liquid  helium  temperature  runs 
was  caused  by  the  plugging  of  the  lower  heat  exchanger  with  solid  hydrogen  - the  result  of 
100  ppm  in  the  helium  gas  supply.  Replacement  of  the  gas  supply  eliminated  this  problem. 

In  interpreting  the  experimental  results  and  comparing  them  to  the  analysis,  we 
should  keep  in  mind  the  strong  temperature  dependence  of  the  test  section  heat  capacity, 
figure  1.11.  This  change  in  heat  capacity,  coupled  with  the  variation  of  the  helium 

4 

density,  changes  b (the  heat  capacity  ratio)  by  a factor  of  about  2(10)  between  300  and 
4 K.  In  spite  of  these  strongly  temperature-dependent  properties,  however,  the  agreement 
of  the  experiment  with  the  analysis  is  quite  good  when  average  property  values  are  used. 

The  response  of  the  wall  temperature  to  a near  step  change  in  coolant  inlet  tempera- 
ture is  shown  in  figures  1.12  through  1.14.  The  slight  asymmetry  in  figures  1.12  and  1.13 
is  due  to  a larger  time  constant  for  the  return  stream  inlet  temperature  compared  to  that 
for  the  go  stream  inlet  temperature.  Single-ended  cool-down  (figure  1.14)  gives,  as  ex- 
pected, temperature  profiles  similar  to  those  for  a half-length  of  the  double-ended  case. 

For  cool-down  to  helium  temperatures  (figure  1.13  and  1.14)  the  curvature  near  the 
cool  ends  is  reduced,  or  becomes  slightly  concave,  compared  to  figure  1.12,  because  of  the 
strong  temperature  dependence  of  the  wall  heat  capacity  at  low  temperature. 

By  using  the  normalized  wall  temperature  ratio  w/w^,  we  may  compare  the  analytical 
(equation  1.14),  numerical,  and  experimental  wall  temperature  profiles,  figure  1.15  during 
the  second  stage  of  cool-down.  The  experimental  results  are  for  cool-down  to  77  K,  so  the 
wall  heat  capacity  varies  by  a factor  of  about  2 in  the  experiment,  whereas  the  numerical 
and  analytical  results  assume  a constant  heat  capacity.  For  w^  < .75,  the  three  results 
essentially  coincide,  and  we  may  expect  the  sinusoidal  approximation  to  accurately  describe 
the  temperature  profile  during  most  of  the  cool-down  if  the  temperature  dependence  of  the 
heat  capacity  is  not  too  large. 

As  discussed  in  Section  1.2,  the  key  to  the  analysis  is  to  obtain  the  near-end  tem- 
perature difference  between  the  go  and  return  streams,  v^  - u^,  since  this  temperature 
difference  determines  the  rate  of  refrigeration,  and  hence,  the  cool-down  rate.  We  may 
compare  the  analytical,  numerical,  and  experimental  values  of  v^  - u^  by  noting  that 
u^  = 0 and  rearranging  1.17  to  give: 


V 

0 


u 


L 


£ 


7t 


1.35 


Using  this  parameter,  we  may  compare  values  of  v^  for  different  values  of  L and  u^. 

Figure  1.15  shows  that  equation  (1.35)  gives  a reasonable  approximation  of  both  the 
numerical  and  experimental  results.  Because  the  experimental  values  of  v^  fall  about  five 
percent  below  equations  (1.17)  and  (1.35),  the  experimental  cool-down  times  should  be 
expected  to  be  about  five  percent  greater  than  predicted  by  the  sinusoidal  approximation. 

Experimental  cooling  curves  in  dimensional  time  and  temperature  are  given  in  figures 
1.16  through  1.18.  The  temperature  is  the  maximum  temperature  in  the  line  (position  x/£  = 1 
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Table  1. --Summary  of  Experiments 


EXP.  NO. 

TEMPERATURE 

(K) 

FLOW 

(g/s) 

COMMENTS 

Counterflow  Steady  State 

1/29 

180 

- 

76 

.2 

N,  = 16.1 
tu 

2/12 

274 

- 

77 

.2 

N.  = 16.1 
tu 

2/12 

257 

- 

77 

.1 

N.  = 21.9 
tu 

3/16 

277 

- 

87 

.3 

N.  = 13.0 
tu 

4/9 

280 

- 

4.2 

.2 

temperature  control  problem 

4/20 

270 

- 

4 

.2 

temperature  control  problem 

4/25 

279 

- 

78 

.2 

N.  = 15.9 
tu 

5/1 

275 

- 

4 

.2 

N.  = 18.2 
tu 

6/27 

277 

- 

76 

.45 

N = 10.0 
tu 

Single  ! 

Stream-Cool  Down 

1/26 

180 

- 

79 

.2 

(pAC)^^  = 1.514  (10)^  J/K 

2/13 

251 

- 

83 

.2 

(pAC)^^  = 1.657  (10)^ 

3/2 

265 

- 

79 

.3 

(pAC)^^  = 1.647  (10)^ 

3/16 

271 

- 

15 

.2 

(pAC)^^  = 1.250  (10)^ 

Counter 

Flow-Cool  Down 

- 

Both 

Ends  Cooled 

1/26 

179 

- 

76 

.2 

2/9 

274 

- 

76 

.2 

2/14 

255 

- 

76 

.1 

3/1 

278 

- 

76 

.3 

3/7 

278 

- 

6 

.2 

aborted  helium  run-H2  contamination 

3/14 

278 

- 

6 

.2 

H II  II  II 

4/3 

274 

- 

7 

.2 

6/28 

281 

- 

75 

.45 

Counter 

Flow-Cool  Down 

- 

One  End  Cooled 

5/8 

270 

- 

6 

.2 

5/9 

2&9 

- 

77 

.2 

5/14 

273 

- 

76 

.3 

5/29 

278 

- 

75 

.45 

Counter 

Flow  Warm-Up  - 

One  End 

Warmed 

5/8 

22 

- 

277 

.2 

5/15 

118 

- 

273 

.2 

5/29 

82 

- 

278 

.45 
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Figure  1.11  Specific  heat  of  stainless  steel. 


Figure  1.12  Test  section  temperature  profile  during  cool-down  to  77  K with  both  ends 
cooled,  experimental  results:  p = 10  bar,  lii  = 0.2  g/s , = 16.1,  UP£  = 16.6  W/s, 

(pAC)^  = 1657  J/K.  The  curves  are  to  guide  the  eye  through  the  experimental  data. 
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Figure  1.13  Test  section  temperature  profile  during  cool-down  to  7 K with  both  ends 
cooled,  experimental  results:  p = 10  bar,  lii  = 0.2  g/s,  = 18.2,  UP£  = 19.3  w/s, 
(pAC)^  = 1250  J/K.  The  curves  are  to  guide  the  eye  through  the  experimental  data. 


Figure  1.14  Test  section  temperature  profile  during  cool-down  to  6 K with  one  end 
cooled,  experimental  results:  p = 10  bar,  lii  = .2  g/s,  = 18.2,  UP£  = 19.3  w/s, 
(pAC)^  = 1250  J/K.  The  curves  are  to  guide  the  eye  through  the  experimental  data. 
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Figure  1.15  Normalized  temperatue  profile  - comparison  of  experimental,  analytical  and 
numerical  (constant  property)  results. 
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SINUSOIDAL  APPROXIMATION. 
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^ A 

- 

□ 

NUMERICAL,  N.  = 50 
tu 

- 
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1 1 

EXPERIMENTAL,  N^^  =16.1 

1 1 

0 
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CO 

Figure  1.16  Dimensionless  inlet  temperature  difference  - comparison  of  experimental, 
analytical  and  numerical  (constant  property)  results  (double  ended  cool-down). 
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t (min) 

Figure  1.17  Experimental  cooling  curves  for 
cooldown  to  76  K - both  ends  cooled.  The 
temperature  is  at  the  mid-point  (x/£  = 0.5) 
of  the  test  section. 


Figure  1.19  Experimental  cooling  curves  for 
cooldown  to  6 K.  The  temperature  is  at  the 
mid-point  (x/£  = 0.5)  of  the  test  section  with 
both  ends  cooled  and  at  the  far  end  (x/£  = 

1.0)  with  one  end  cooled.  For  both  cases 
rti  = 0.2  g/s. 


Figure  1.18  Experimental  cooling  curves  for 
cooldown  to  76  K - one  end  cooled.  The 
temperature  is  at  the  far  end  {x/l  = 1.0)  of 
the  test  section. 


Figure  1.20  Experimental  warming  curves  for 
warm-up  from  22,  82,  and  118  K with  one  end 
warmed.  The  temperature  is  at  the  far  end 
(x/£  = 1.0)  of  the  test  section. 
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for  one  end  cooled  or  position  x/jZ  =0.5  for  both  ends  cooled).  The  curves  for  cool-down 

to  helium  temperature  are  more  nearly  linear  than  exponential  because  of  the  rapid  decay 

in  the  wall  heat  capacity  at  low  temperatures. 

Experimental  warm-up  curves  are  given  in  figure  1.20  for  single-ended  warming.  For 

constant  properties,  warm-up  is  simply  the  inverse  of  cool-down. 

2 

The  analysis  suggests  the  use  of  x'/N^^  as  a dimensionless  correlating  parameter  for 

the  cool-down  curves.  Then  according  to  equation  (1.26)  or  (1.30)  the  normalized  wall 

1 4 

temperature  w will  exhibit  only  a weak  dependence  ([1  + ] for  doubled  ended  cooling)  on 

^tu 

^tu  Figure  1.21  compares  the  analytical  and  numerical  results  for  zero  heat 

leak  in  this  fashion.  When  adjusted  for  b and  ^ a>,  the  numerical  and  analytical  results 
agree  to  better  than  one  percent.  Considering  the  approximate  nature  of  the  analysis,  this 
is  probably  fortuitous. 

2 

The  dependence  of  the  cool  down  parameter  x'/N^^  on  1/N^^  is  shown  in  figure  1.22,  and 

it  appears  that  the  sinusoidal  approximation  underestimates  this  dependence.  For  N.  >5, 
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t'/Ntu  (1  + j;i^ — ) more  closely  describes  the  numerical  results  (both  ends  cooled). 

2 

The  dependence  of  x'/N^^^  on  the  heat  capacity  ratio  b is  shown  in  figure  1.23.  In  this 
case  the  analysis,  x'/N^^  (1  + ^),  agrees  with  the  numerical  results  within  the  accuracy 

of  the  calculations.  This  should  be  expected  since  (1  + I)  is  the  total  dimensionless  heat 

D 2 

capacity  for  the  system  described  by  equations  (1.4)  through  (1.6),  i.e.,  (1  + ^)  is  the 

dimensionless  sum  of  the  wall  heat  capacity  plus  the  go  and  return  stream  heat  capacities. 

Thus,  for  this  simple  system,  the  cool-down  time  is  proportional  to  the  total  heat  capacity 

of  the  walls  and  fluid,  but  independent  of  its  distribution. 

On  the  basis  of  figures  1.21  and  1.22,  we  may  modify  the  analytical  results  to  give  a 

x'/N^  (1  + ~^  ) dependence.  Equations  (1.26),  (1.29),  and  (1.30)  then  become: 

^tu 


~1  - 1 

exp 

7t  t2 

, QL  1.26a 

7t  me  (T.  - T.) 

L p 1 f 0 

2L^  (1  + nh  (1  + h 

L 4L  ^ J 

a *Cp(T.  - ip 

X',  = 1.048  (1*^)  (l.|)  1.29a 

for  single  ended  cool  down,  and 

for  doubled-ended  cool-down  (90  percent  cool  down  for  both  cases).  Equations  (1.26a), 

1.29a),  and  (1.30a)  represent  the  numerical  constant  property  and  mass  flow  results  to  within 
one  percent  for  N^^IO  and  to  within  about  four  percent  for  5<N^^<10. 
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1/N 


tu 


Figure  1.21  Analytical  and  numerical  (constant  property)  cooling 
curves  in  terms  of  of  dimensionless  parameters. 


Figure  1.22  Dependence  of  the  cool-down  parameter  on  1/N^,,.  The  parameter 


2 

x'/Ntu  is  evaluated  at  w = 0.05. 


tu‘ 
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2 

Figure  1.23  Dependence  of  the  cool-down  parameter  I'/N,  on  the  heat  capacity  ratio 

2 « u 

b.  The  parameter  is  evaluated  at  w = 0.05. 


Figure  1.24  Comparison  of  the  experimental  and  modified  analytical  (equation  1.22a) 
cooling  curves  for  cool  down  to  76  K.  The  results  are  given  in  dimensionless 
parameters. 
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Figure  1.25  Experimental  cooling  curves  in  dimensionless  parameters  for  cool-down  to 
6 K.  The  normalized  temperature  w is  given  at  x/£  =0.5  for  both  ends  cooled  and  at 
x/£  =1.0  for  one  end  cooled. 


Figure  1.26  Experimental  warming  curves  in  dimensionless  parameters. 
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A comparison  of  (equation  1.26a)  with  the  LN^  temperature  cooling  experiments  is 

^ 2 

given  in  figure  1.24.  Note  that  the  parameter  coalesces  the  experimental  data  to  a 

single  curve,  even  though  the  dimensional  cool-down  times  (figure  1.17  and  1.18)  differ  by 
more  than  a factor  of  10.  To  generate  the  curve,  equation  (1.26a)  was  evaluated  for 
= 15  and  QL/itCp  (T^.  - T^)  = 0.05,  values  typical  of  the  experiments.  There  are  several 
differences  between  experiment  and  analysis.  In  the  experiment,  the  test  section  heat 
capacity  varies  with  temperature.  The  heat  leak  varies  both  with  time  during  an  individual 
experiment  and  varies  from  experiment  to  experiment  as  well.  Considering  these  differences 
between  experiment  and  analysis,  the  agreement  of  the  experiment  with  the  analysis  is 
better  than  expected. 

Figure  1.25  presents  the  helium  temperature  cooling  curves  (figure  1.18)  as  a function 
of  the  cool-down  parameter  As  with  the  liquid  nitrogen  temperature  experiments, 

this  formulation  coalesces  the  data  to  a single  curve.  The  convex  shape  of  the  curve  is 
due  to  the  large  variation  in  the  test  section  heat  capacity  with  temperature.  Complete 
cool  down  is  achieved  at  ^ 0.3  (double  ended  cooling). 

Comparison  of  the  warm  up  experiments  with  the  cool-down  experiments  (figure  1.26) 
shows,  as  expected,  that  warm  up  is  essentially  the  inverse  of  cool-down.  Initially, 
warm-up  proceeds  somewhat  faster  because  of  the  lower  initial  heat  capacity.  The  effect 
of  heat  leak  is  much  larger  in  the  warm-up  experiments  since  the  final  stage  of  the  warm- 
up occurs  at  ambient  temperature  where  radiation  heat  transfer  effects  are  much  larger. 

In  the  warm-up  of  a real  SPTL,  the  heat  leak  will  always  be  positive  (into  the  line)  so 
that  heat  leak  will  aid  the  warm-up.  In  these  warm-up  experiments,  the  heat  leak  is 
negative  and  consequently  slows  down  the  warm-up. 


1.5  Detailed  Computer  Program  for  the  Calculation  of  Cool-Down  - Comparison  with  Experi- 
mental Results 

Previous  to  the  work  reported  in  the  sections  thru  1.4,  one  of  us  (MCJ)  developed  a 
more  detailed  computer  program  for  the  calculation  of  flow  and  heat  transfer  transients  in 
channels  in  which  single  phase  (supercritical)  helium  is  a coolant.  This  program  differs 
from  the  numerical  solution  of  equations  1.4  thru  1.7  in  that  it  considers  energy,  momentum, 
and  continuity  for  each  channel  rather  than  just  energy,  and  in  that  variable  thermal 
properties  are  considered  for  both  the  cable  and  the  coolant.  The  main  features  of  the 
program  are; 

i)  Time-dependent  fluid  behavior  is  described  by  a transformed  set  of  three  partial 
differential  equations  of  conservation  with  dependent  variables  pressure  P,  enthalphy  H, 
and  velocity  U.  These  variables  represent  cross-sectional  averages,  the  flow  being  treated 
one-dimensional ly. 

ii)  Fluid  properties  and  other  thermodynamic  variables  such  as  temperature  are 
represented  by  polynomials  with  P and  H as  arguments.  The  use  of  polynomials  permits  fast 
computation  with  acceptable  error.  These  polynomials  are  contained  in  a properties  sub- 
routine for  continuous  up-dating  to  current  local  conditions. 


25 


iii)  Linear  heat  transfer  to  the  channel  wall  is  allowed  and  the  temperature  of  the 
channel  wall  is  described  by  the  transient  heat  conduction  equation  for  one-dimensional 
heat  flow. 

iv)  The  resulting  set  of  8 partial  differential  equations  is  solved  with  appropriate 
boundary  and  initial  conditions  by  an  available  software  package  based  on  finite  element 
col  location. 

Further  details  can  be  found  in  reference  [1]. 

The  computer  program  has  been  used  to  make  predictive  calculations  on  the  cool-down 
of  superconducting  power  transmission  lines  (SPTL's)  and  for  the  study  of  transients  in 
hollow  superconductors.  Now  that  experimental  results  are  available  for  cool-down,  it 
seemed  wise  to  validate  the  computer  program  by  comparison  of  computed  with  experimental 
results.  This  work  is  now  partially  complete  and  here  we  report  the  results  of  a compari- 
son with  an  experiment  in  which  a conduit  with  two  thermally  coupled  streams  in  counter 
flow  was  cooled  with  supercritical  helium  to  liquid  nitrogen  temperature.  The  discussion 
following  refers  to  the  run  of  3/1/79  in  which  a mass  flow  rate  of  0.3  g/s  was  fed  in  at 
each  end  of  the  test  section. 

The  comparison  we  have  made  shows  first  of  all  that  perhaps  the  least  well-known 
quantity  is  the  time  dependent  heat  leak.  The  experiments  described  in  this  report  (see 
figure  1.10)  show  that  this  is  a complex  function  since  it  first  increases  and  then  de- 
creases with  time  due  to  experimental  conditions  discussed  in  Section  1.4.  It  was  there- 
fore necessary  to  fit  the  heat  leak  function  to  the  experimental  results.  It  was  found 
that  for  all  runs  a single  crude  approximation  of  the  form 

Q,eak  = 1.36 

could  be  used,  where  Qig^j^  is  a lumped  heat  leak  and  a,b,c,  and  d are  empirical  constants. 

Results  for  the  mid-point  temperature  computed  as  a function  of  time  using  the  above 
heat  leak  function  are  shown  in  figure  1.27.  In  the  dashed  curve,  the  actual  measured 
mass  of  the  two-stream  channel  was  used.  In  all  other  respects  this  was  a calculation 
from  first  principles  and  while  the  error  is  acceptable  for  engineering  purposes  it  indi- 
cates an  imprecise  knowledge  of  the  thermal  capacity  of  the  two-stream  channel  or  the  heat 
transfer  resistance  between  the  two  streams,  or  both.  Now  the  actual  thermal  capacity  is 
known  experimentally  from  single-stream  cool-down  experiments.  In  that  determination  (and 
in  the  cool-down  computation)  the  heat  capacity  data  from  [11]  was  used  and  it  was  then 
found  that  the  effective  mass  of  the  test  section  was  4.603  kg.  In  separate  steady-state, 
it  was  determined  that  the  product  of  overall  heat  transfer  coefficient  and  heat  transfer 
area  (DA)  was  20.31  W/K.  These  figures  compare  with  the  measured  mass  of  4.095  kg  and  a 
calculated  UA  of  19.54  W/K  used  in  the  computation.  We  therefore  ran  the  computation 
again  with  the  larger  figures.  The  result  for  mid-point  temperatures  is  the  solid  curve 
in  figure  1.27.  The  agreement  with  experiment  is  now  seen  to  be  almost  exact. 

A comparison  of  computed  and  measured  temperature  profiles  at  two  different  times 
during  cool-down  is  shown  in  figure  1.28.  Again  the  agreement  is  almost  exact. 
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Figure  1.27  Comparison  of  computed  mid-point  temperature  with  experimental  cool-down 
results.  0,  experimental; — computed  (measured  test  section  mass); computed  (effec- 

tive test  section  mass). 


Figure  1.28  Comparison  of  computed  with  experimental  temperature  profiles.  Experimental 
results,  • 960s,  o 4560s.  Computed  curves  are  for  1000s  and  4500s  respectively. 
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With  the  comparison  at  SPTL  operating  temperature  yet  to  be  completed,  we  conclude 
tentatively  that  the  computer  program  is  capable  of  excellent  prediction  provided  the  heat 
leak  function,  the  effective  mass,  and  the  heat  transfer  resistance  are  accurately  known. 

The  latter  two  quantities  can  apparently  be  estimated  from  first  principles  - at  least  for 
the  simple  experimental  simulation  - with  acceptable  accuracy.  Because  of  its  relatively 
high  value  in  the  laboratory  experiments,  the  heat  leak  function  is  more  critical  and  in 
this  instance,  it  was  necessary  to  rely  on  the  experimentally  determined  function.  In  the 
case  of  cool-down  calculations  for  full-scale  SPTL's,  the  relative  size  of  the  heat  leak 
is  much  smaller  and  the  results  are  not  particularly  sensitive  to  the  heat  leak  function. 

1.6  Summary 

We  have  investigated  the  two  stream  counterflow  method  of  SPTL  cool-down  by  a com- 
bined analytical,  numerical,  and  experimental  approach.  The  experiments  were  conducted 
with  wel 1 -control  1 ed  boundary  conditions  and  constant  mass  flow  rate  for  ease  of  com- 
parison with  the  analytical  and  numerical  methods. 

Counter  flow  cool-down  progresses  in  two  stages.  During  the  first,  and  shorter, 
stage  a constant  amplitude  temperature  wave  develops  until  it  extends  the  full  length  of 
the  line,  figure  1.3.  During  the  second  stage  the  amplitude  of  this  wave  decays  exponen- 
tially to  zero,  figure  1.4.  The  rate  of  decay  is  proportional  to  the  temperature  difference 
between  the  go  and  return  streams  (i.e.,  to  the  rate  of  refrigeration)  which  in  turn  is 
proportional  to  the  thermal  resistance  between  the  go  and  return  streams,  i.e. 

I ^ U(-)^ 

for  large  Thus  for  perfect  heat  transfer  between  the  two  streams  (U->a>)  the  cool-down 

time  goes  to  infinity.  Clearly,  poor  thermal  contact  between  the  go  and  return  streams  is 
required  for  short  cool-down  times. 

The  flow  rate  in  a full-scale  transmission  line  will  usually  be  limited  by  the  avail- 
able pressure  drop,  in  which  case  analysis  gives 

cd  pAp  . 


The  more  detailed  results  of  the  analysis  are  given  by  equations  1.26a,  1.29a  and 
1.30a. 

We  conclude  by  noting  that  the  analytical  results  give  a reasonable  prediction  of  the 
cool-down  time  for  the  constant  mass  flow  rate  case  when  average  property  values  are  used. 
The  detailed  numerical  model  reproduces  the  experimental  results  for  cool-down  to  80  K so 
that  we  may  expect  that  it  will  also  give  accurate  estimates  for  cool-down  to  SPTL  operat- 
ing conditions  for  a variety  of  boundary  conditions. 
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1.7 


Nomenclature 


A 

b 


C 

C 

D 

f 

h 

£ 

L 

rti 

N 

P 

P 


P 


tu 


Q 


t 

T 

u 

U 

V 

V 
w 

X 

y 


cross  sectional  flow  area 
heat  capacity  ratio  = 

specific  heat  of  solid 

specific  heat  at  constant  pressure  for  gas 
hydraulic  diameter 
friction  factor 

film  coefficient  + 1/2  solid  resistance 
length  of  transmission  line 
dimensionless  length  = 

mass  flow  rate  of  gas 

u .c  ^ .c  -a.  UP£ 
number  of  transfer  units  = xtt- 

liiCp 

pressure 

heated  perimeter 

total  heat  leak  for  single  ended  cool  down, 

^2  the  total  heat  leak  for  double  ended  cool  down. 

time 

temperature 

normalized  temperature  of  go  stream  = (T  - T^)|^/(T^.  - 

overall  heat  transfer  coefficient  between  the  go  and  return  streams. 

normalized  temperature  of  return  stream  = (T  - Tp)^/(T^.  - 

gas  or  fluid  velocity 

normalized  wall  temperature  = (T  - T^)^/(T.  - 
position 

h Px  X 

dimensionless  position  = tjt-  = L-r 

^ rtiCp  £ 


Greek 

e heat  exchanger  effectiveness 

dimensionless  position  of  wave  front  during  the  first  phase  of  cool  down. 


density  (gas  if  unsubscripted) 
hP  . 


dimensionless  time  = 
dimensionless  time  = 


(pAC) 

UP 

(pAC) 


Subscripts 
cd  cool  down 

i initial  condition 

in  inlet 

f final  condition  for 

£ X = £ position 

0 X = 0 position 


out 

u 

V 

zero  heat  input  w 

1 

2 


outlet 
go  stream 
return  stream 
wal  1 

first  phase  of  cool  down 
second  phase  of  cool  down 
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1.9  Appendix  - Cool-down  time  for  constant  pressure  drop  boundary  condition 


Starting  with  the  general  expression  for  pressure  drop  in  a tube, 


dp  = - 


32f 


P 


dx , 


and  assuming  an  ideal  gas  equation  of  state, 


P = E- 


RT 


gi  ves 


dp  = - 


32fR 


where  T = T(x)  and  rii  = m(t). 

The  sinusoidal  temperature  distribution  assumed  in  equation  1.13  may  also  be 
as 


T(x)  = (T.  - T ) sin  n x + T 
£0  2£  ° 


Substituting  1.3a  into  1.2a  and  integrating  over  the  length  of  the  channel  gives 


lii^  = pAp 


£ C[2  (T  - T ) + T ] 

- £ o 0^ 

n 


where 


P + P 
p = £ o 


AP  = P - P„ 
o £ 


C = 32fR 


71^0" 


To  obatin  the  cool-down  time  for  constant  pressure  drop  in  the  go-stream,  we 
to  equation  1.12,  the  first  law  of  thermodynamics  for  the  transmission  line; 


/' 


mC  v^dt  = - (pAC)^(l  + 2)  /[w(x,t  + dt)  - w(x,t)]dx  . 

p o b 


1.1a 


1.2a 


wri tten 


1.  3a 


1.4a 


return 


1.5a 
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The  normalized  temperature  difference  between  the  go  and  return  streams,  v^,  is  given 


by  equation  1.17  as 


or,  noting  that 


V = 7t  u„ 
o L ^ 


L = hP£ 
liiCp 


V 

o 


= nrtiCp 
hP 


1.6a 


Substituting  1.6a  into  1.5a,  replacing  m^  with  1.4a  and  evaluating  the  integral  on  the 
right  hand  side  of  1.5a  gives 


n2Cp2  PAP  dt  = - du.(pAC)  (1  + 2)(1  + n) 

hP  ' - T ) + T ] b 2L 


Noting  that 


T ) + T 
o 0 


-^£ 

71 


equation  1.7a  becomes 


/ 2 + a \ du,  = - B PAP  dt 


where 


a 


T 

o 


1.7a 


1.8a 


B = Ti^Cp^D^  1.9a 

64(pAC)^(l  + 2)(1  + 7t)hPfR(T.  - T^) 

w b 2L  10 

If  we  neglect  the  time  dependence  of  L,  ie  assume  that  L is  large,  then  we  may  integrate 
1.9a  to  give 

2(u^  + 1)  +alog(Uj,l"  B / PAP  \ t 1.10a 


for  the  initial  condition  of  u^  = 1 at  t = 0. 
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2.0  CABLE  PERMEABILITY 
(P.  R.  Ludtke  and  D.  E.  Daney) 


In  case  of  a rupture  in  the  lead  sheath  of  the  cable,  it  is  important  to  know  how  much 
leakage  from  the  pressurized  (15  bar)  cable  might  occur.  Our  measurements  of  the  cable 
permeability  were  intended  to  define  the  magnitude  of  the  problem  rather  than  to  give  a 
precise  value  for  the  leak  rate  at  the  anticipated  operating  conditions. 

The  apparatus  and  the  cable  are  the  same  ones  that  were  used  previously  for  thermal 
cycling  of  a BNL  model  cable  [1].  .Figure  2.1  is  a cut  back  view  of  the  cable.  As  shown 
in  figure  2.2,  five  symmetrically-spaced  gas  taps  in  the  cable's  lead  sheath  are  connected 
to  individual  wet  test  gas  flow  meters.  The  taps  are  formed  by  soldering  (Wood's  metal 
solder)  4.8  mm  ID  tubes  into  holes  cut  in  the  lead  sheath.  Thus,  the  experiment  measures 
the  leak  rate  through  the  composite  between  the  pressurized  bronze  core  of  the  cable  and 
the  4.8  mm  diameter  holes  in  the  lead  sheath  spaced  at  70  cm  intervals.  The  effect  of  end 
leakage  (between  the  core  and  the  lead  sheath)  is  given  by  the  difference  between  the  flow 
rates  of  the  center  tap  and  the  end  taps.  The  ambient  pressure  is  .82  bar,  and  the  pres- 
sure difference  across  the  cable  is  the  gauge  pressure.  Helium  gas  was  used  in  all  tests. 

Figures  2.3  through  2.5  give  the  measured  cable  leak  rate  at  ambient  temperature  and 
at  76  K.  Some  shift  in  the  ambient  leak  rate  occurred  even  though  the  cable  had  previously 
been  thermally  cycled  20  times. 

In  order  to  estimate  the  leak  rate  at  operating  conditions  (P  = 15  bar  and  T = 8 K), 
we  use  the  expression: 


and  obtain  a leak  rate  of  about  15  g/s  per  hole  - extremly  high.  Here  G is  mass  flux, 
p density,  AP  pressure  difference,  and  p viscosity. 

From  these  measurements  we  conclude  that  any  significant  rupture  of  the  lead  sheath 


ability  of  the  composite  within  the  lead  sheath  also  suggests  that  some  fluid  convection 
in  this  layer  may  accompany  forced  flow  inside  the  bronze  core  since  the  composite  layer 
will  experience  the  same  axial  pressure  gradient.*  As  a result,  the  thermal  conductance 
across  the  wrapped  dielectric  may  be  very  much  higher  than  anticipated.  Because  this 
conductance  is  critical  to  the  performance  of  the  SPTL,  we  strongly  recommend  that  the 
radial  thermal  conductance  of  a cable  be  measured  under  the  expected  operated  conditions, 
i.e.,  with  supercritical  helium  flowing  through  the  core. 


*Natural  convection  effects  in  wrapped  dielectrics  have  been  measured  separately  by 
Thullen  and  Hajjar  (F.  G.  Hajjar,  Thermal  Conductivity  of  Selected  Dielectric  Materials  in 
Liquid  Nitrogen,  paper  AIAA-80-0089 , AIAA  18th  Aerospace  sciences  meeting,  January  14-16, 
1980,  Pasadena,  California). 


will  lead  to  an  unacceptably  high  leak  rate  from  the  pressurized  core.  The  high  perme- 
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Figure  2.1  Cable  cross  section,  cut  back  view 
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6.35  mm  O.D. 


Figure  2.2  Cable  permeability  test  arrangement. 


GAUGE  PRESSURE,  Pa 


Figure  2.3  Cable  leak  rate  before  cool-down. 


35 


FLOW 


Figure  2.4  Cable  leak  rate  at  76  K. 


Figure  2.5  Cable  leak  rate  after  cool-down. 
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3.0  EVALUATION  OF  THERMAL  FLUX  METERS 
(P.  R.  Ludtke) 


3.1  Introduction 

Integrity  of  the  insulating  vacuum  in  the  transmission  line  thermal  enclosure  is 
essential  to  the  successful  operati'on  of  the  line.  The  standard  practice  for  determining 
the  integrity  of  a multilayer  insulated  cryostat  is  to  measure  the  insulating  vacuum 
directly.  The  required  penetration  of  the  vacuum  space  as  well  as  the  vacuum  gauge  itself, 
however,  introduce  an  additional  risk  of  a vacuum  leak.  For  this  reason,  it  is  highly 
desirable  to  determine  the  integrity  of  the  thermal  insulation  by  some  non-intrusive 
means. 

Thermal  flux  meters,  which  are  the  subject  of  this  study,  have  been  proposed  for  this 
purpose.  Not  only  are  these  meters  non-intrusive  because  they  are  mounted  on  the  outside 
of  the  enclosure  wall,  but  they  also  provide  the  enclosure  heat  flux  directly. 

Successful  field  application  of  heat  flux  meters  requires  detailed  study  of  a number 

of  factors  such  as  meter  sensitivity  and  effects  of  moisture  and  other  soil  conditions. 

Most  important,  however,  are  the  variations  in  the  soil  heat  flux  which  might  mask  the 

2 

rather  low  enclosure  heat  leak  - approximately  30  pW/cm  . These  variations  in  the  soil 
heat  flux  have  three  basic  periods;  the  diurnal  variation,  variations  of  a few  days  or 
weeks  due  to  weather  patterns,  and  the  seasonal  (yearly)  variation.  Since  the  amplitude 
of  the  variation  in  soil  temperature  is  diminished  by  the  factor 

exp  (-  ^n/ax^  x) 

the  effect  of  the  seasonal  variation  will  be  far  larger  than  the  diurnal  variation.  In 
the  above,  a,  is  the  thermal  diffusivity  of  the  soil,  is  the  period  and  x is  the  depth. 

Preliminary  concern  was  centered  on  the  diurnal 
variation,  and  analysis  indicated  a signal  to  noise  ratio  of  over  5,000  if  the  sensors  are 
buried  at  a depth  of  1.5'meters.  The  results  so  far  indicate  that  the  effect  of  variable 
soil  heat  flux  can  largely  be  eliminated  by  taking  the  algebraic  sum  of  the  signals  mounted 
180°  from  each  other. 

During  this  reporting  period,  laboratory  testing  of  the  heat  sensors  was  completed, 
and  field  testing  was  initiated. 


3.2  Experimental 

The  experimental  study  proceeded  in  two  stages.  First,  the  sensitivity  of  four  pairs 
of  commercial  heat  flux  sensors  were  evaluated  in  the  laboratory.  Then,  the  sensors  were 
evaluated  under  field  conditions,  i.e.,  buried  and  subjected  to  diurnal  and  longer  term 
variations  of  the  soil  temperature  and  heat  flux. 
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The  four  pairs  of  heat  sensors  from  different  manufacturers  are  mounted  on  a 51  cm 
(20-inch)  O.D.  x 61  cm  (24-inch)  long  steel  pipe  which  simulates  the  cryogenic  envelope. 
Table  3.1  summarizes  the  properties  of  the  meters.  Six  of  the  meters  have  flat  surfaces, 
and  the  other  two  have  surfaces  contoured  to  fit  the  outer  surface  of  the  pipe.  Copper 
adapters  with  one  side  flat  and  the  other  side  curved  were  machined  from  high  purity 
copper.  The  adapters  were  designed  such  that  the  curved  surface  area  was  equal  to  the 
flat  surface  area  of  the  heat  flow  meters.  A moderately  high-thermal -conduct! vi ty  epoxy 
was  used  to  bond  the  adapters  and  meters  to  the  steel  pipe.  Sensors  from  each  pair  were 
mounted  180°  from  each  other  on  a circumferential  centerline  equidistant  from  each  end  of 
the  pipe.  The  outer  surface  of  the  pipe  was  divided  by  four  quadrant  lines  parallel  to 
the  pipe's  axis.  Two  sensors  were  placed  near  each  quadrant  line;  the  edge  of  each  sensor 
was  2 cm  from  the  quadrant  line. 

Early  in  the  program,  four  sandwich-type  meters  delaminated.  The  copper  or  aluminum 
surface  plates,  which  are  bonded  to  the  inner  plastic  core,  separated  and  three  of  the 
four  meters  were  replaced  by  manufacturers  C & D prior  to  testing.  Later  in  the  program, 
the  fourth  sensor  malfunctioned  28  days  into  the  field  test.  Meters  A and  B are  of  one 
piece  plastic  construction,  and  appear  to  be  much  better  sealed  for  an  underground  ap- 
plication where  moisture  is  of  concern.  All  of  the  meters  have  a monitoring  thermocouple 
within  the  meter  body. 

Uniform  heat  flux  is  generated  by  a heater  mounted  on  the  inside  surface  of  the  pipe. 
It  is  fabricated  from  60  series-connected  segments  of  stainless  steel  shim  stock  .05  mm 
(.002-inch)  x 25.4  mm  (1-inch)  x 60  cm  (23  3.4-inch),  bonded  by  RTV  silicone  to  .13  mm 
(.005-inch)  polyimide  film  electrical  insulation. 


Table  3. 1--Heat-Fl ow  Meter  Specifications 


Si  ze 

Contoured 

Manufacturer' s 

Type  of 

Type  of 

Brand 

(cm) 

Surface 

Sensi ti vi ty 

Construct! on 

Sensor 

watts 

- m\/ 

A 

11.4 

X 

11.4 

X .48 

NO 

300 

Integral , 
Epoxy  Glass 

Thermopi 1 e 

B 

5.08 

X 

15.2 

X .64 

YES 

135 

Integral , 
Diallyl 
Phthal ate 

Thermopi 1 e 

C 

15.2 

X 

15.2 

X .47 

NO 

113 

Sandwi ch , 

Solid  State 

Copper- 

Polyimide 

Thermopi 1 e 

D 

5.08 

X 

5.08 

X .25 

NO 

160 

Sandwi ch , 

Solid  State 

A1 umi num- 
Phenol ic 

Thermopi 1 e 
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Potential  taps  at  the  heater  ends  are  used  for  measuring  the  electric  power  input. 

The  effects  of  end  losses  are  minimized  by  mounting  the  heat-flux  sensors  at  the  axial 
center  of  the  pipe.  A photo  of  the  heater  strips  inside  the  pipe  section  is  shown  in 
figure  3.1.  The  inside  of  the  pipe  was  then  filled  with  pour-type  polyurethane  foam,  and 
the  ends  were  hermetically  sealed  with  fiberglass  caps.  A thin  coat  of  RTV  silicone 
coating  was  applied  to  the  sensors  with  laminar  construction  for  an  additional  hermetic 
seal.  Photos  of  the  complete  pipe  section  assembly,  with  heat-flow  meters,  adapter 
plates,  and  fiberglass  end  caps  is  shown  in  figure  3.2. 

The  laboratory  tests  were  conducted  in  an  isothermal  environment  provided  by  the 
arrangement  in  figure  3.3.  The  temperature  of  the  copper  thermal  shield  is  controlled  to 
± .01°C  by  water  circulated  from  the  temperature  controlled  bath.  Additional  stability  is 
provided  by  the  outer  layer  of  perlite  and  styrofoam  insulation. 


3.2.1  In-situ  Calibration 

Most  of  the  measurements  were  made  with  air  in  the  annulus  surrounding  the  pipe. 
Because  the  field  tests  were  to  be  made  with  the  pipe  buried  in  sand,  a few  laboratory 
tests  were  also  made  with  the  annulus  filled  with  sand.  Somewhat  to  our  surprise,  we 
observed  moderate  shifts  in  the  calibrations  in  going  from  air  to  sand. 

The  expected  heat  flux  to  a typical  cryogenic  enclosure  was  assumed  to  be  that  for 

the  BNL-LILCO  plan  3B  - about  28  pW/cm^.  Accordingly  the  meters  were  tested  at  9,  18,  28, 

53,  114,  and  969  pW/cm  . After  the  pipe  appeared  to  reach  equilibrium  at  a given  power 

level,  measurements  were  continued  for  4 to  8 days  to  obtain  equilibrium  data.  The  emf 

from  each  heat-flow  meter  is  recorded  daily  during  the  equilibrium  period  and  from  these 

data  points,  the  arithmetic  means,  the  sample  standard  deviation,  and  the  coefficient  of 

variation  are  calculated.  Table  3.2  shows  data  from  a typical  calibration  test,  and 

figure  3.4  is  a plot  of  the  coefficients  of  variation  for  the  first  eight  calibration 

2 

tests,  with  air  surrounding  the  meters.  The  highest  heat  flux  of  969  pwatts/cm  has  an 
extremely  low  coefficient  of  variation,  and  all  the  meters  show  a fairly  consistent  coef- 
ficient of  1 to  3%. 

The  temperature  dependence  of  the  meter  calibration  was  determined  at  a power  level 
of  28  pW/cm  and  temperatures  of  4,  15,  and  26°C--the  temperature  range  anticipated  during 
the  field  test.  The  results  of  these  tests,  which  are  given  in  table  3.3,  indicate  a 
slight  dependence,  the  sign  of  which  varies  with  the  meter  type.  A typical  manufacturer's 
correction  curve  is  given  in  figure  3.5.  It  indicates  that  the  temperature  correction 
factor  is  negligible  for  small  deviations  from  the  calibration  temperature. 

The  results  of  the  air  and  sand  calibration  are  given  in  figures  3.6  through  3.9.  In 
all  cases  the  meter  response  is  quite  linear,  but  the  zero  bias  is  somewhat  larger  with 
sand.  The  difference  between  the  air  and  sand  calibration,  summarized  in  table  3.4,  is 
rather  puzzl ing--both  because  it  occurs  at  all,  and  because  it  is  in  different  directions 
for  different  meters. 

Several  factors  could  be  responsible,  such  as  meter  construction,  adapter  plate 
design,  and  pressure  or  moisture  sensitivity.  Meters  A and  B,  which  exhibit  a decrease  in 
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Figure  3.1--Heater  mounted  on  inside  of  diameter  of  pipe. 


Figure  3.2--Heat  flux  mounted  on  pipe  and  ready  for  testing. 
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Figure  3.3  Heat  sensor  laboratory  test  apparatus. 
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Figure  3.4  Coefficient  of  variation  for  laboratory  evaluation  of  heat  flux  meters. 
Different  symbols  are  for  different  heat  fluxes. 
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Table  3. 2--Cal ibration  Test  No.  3 


(Heat  Flux  = 28.5  pWatts/cm^) 


Data  ^ 24  - 168  Hrs. 


ve 


6-2 

11am 

0 


6-3 

6-4 

6-4 

6-5 

6-5 

6-6 

6-6 

6-7 

6-7 

6-8 

6-9 

11am 

9am 

4pm 

8am 

5pm 

8am 

3pm 

9am 

4pm 

4pm 

11am 

24 

46 

53 

69 

78 

93 

100 

118 

125 

149 

168 

Mean 

(m) 


Standard  Coefficient 

Deviation  of  Variation 
s X 100 


(s) 


m 

(%) 


Heat-Flow  Meter  emf  ((jVolts) 


Ijvolts  jjvolts 


176 

183 

550 

555 

487 

539 

491 

538 


96 
101 
306 
308 
272 
301 
270 
295 
Equi 1 ibrium 


96 

95 

101 

96 

100 

98.5 

100 

100 

101 

101 

100 

106 

100.5 

104 

103 

104 

103 

104 

307 

301 

319 

304 

314 

312.5 

320 

317 

321 

309 

303 

322 

304 

317.5 

314 

320 

316 

321 

275 

270 

282 

270 

277.5 

277 

280 

277 

282. 

303.5 

300 

314 

300 

309 

307 

310.5 

309 

314 

270 

266 

281 

266 

277 

274 

279 

274 

281 

295.5 

289 

307 

290 

302 

300 

306 

303 

306 
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102 

105 

320 

320 

277 

309 

281 

307 


98.68 

102.86 

312.86 

314.05 
276.36 
307.0 
274.45 

300.05 


2.49 

1.98 

7.25 

6.95 

4.32 

5.18 

5.82 

6.70 


2.52 

1.92 

2.32 

2.21 

1.56 

1.69 

2.12 

2.23 


Table  3.3--Meter  Dependence  on  Operating  Temperature 


Meter 

Temperature 

4°C 

15°C 

26°C 

Test  No. 
Heat  Flux 

5 

28.8 

cm'^  - 

3 

28.4 

6 

28.4 

cm^ 

Meter  No. 

Mean  Heat-Flow  Meter  emf  (pV) 

A-1 

97.17 

98.68 

101 

A-2 

102.47 

102.86 

103.6 

B-1 

311.41 

312.86 

314.13 

B-2 

313.69 

314.05 

313.75 

C-1 

283.36 

276.36 

268.7 

C-2 

313.67 

307.00 

298.25 

D-1 

278.23 

274.45 

270.5 

D-2 

302.79 

300.05 

297.0 

Table  3. 4--Heat-Flow  Meter  Sensitivities 

Sensor 

Manufacturer's  In-Situ 

Sensitivity  (Air) 

In-Situ 

(Sand) 

A-1 

328.3 

cm^-mV 

323.5 

cm'^-mV 

371.5 

cm'^-mV 

A-2 

287.3  " 

307.7  " 

361.7  " 

B-1 

131.5  " 

98. 

105.2  " 

B-2 

137.3  " 

100.4  " 

110.8  " 

C-1 

117.6  " 

111.5  " 

108.5  " 

C-2 

108.9  " 

99.6  " 

79.8  " 

D-1 

167.3  " 

112.8  " 

99.8  " 

D-2 

154.8  " 

102.4  " 

92.8  " 
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MULTIPLYING  COEFFICIENT 


Figure  3.5  Typical  heat  flux  meter  temperature  correction  curve. 
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Figure  3.6  Calibration  for  sensors  A1  and  A2. 
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SENSOR  emf,  yV 


Figure  3.7  Calibration  for  sensors  B1  and  B2. 


Figure  3.8  Calibration  for  sensors  Cl  and  C2. 
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SENSOR  emf , yV 


Figure  3.9  Calibration  for  sensors  D1  and  D2. 
T ype  E 


thermocouples 
in  sand 


All  dimensions  in  cm. 

Figure  3.10  Burial  configuration  for  the  pipe-meter  assembly. 
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sensitivity,  are  one  piece  cast  plastic,  whereas  meters  C and  D,  which  exhibit  an  increase 
in  sensitivity,  are  of  sandwich  construction  with  metal  plates  for  outer  surfaces.  Both  C 
and  D have  solid-state  sensor  elements.  Heat  loss  through  the  edges  of  the  adaptor  plates 
of  meters  A,  C,  and  D would  cause  a sensitivity  shift  only  if  the  relative  loss  changes 
from  air  to  sand. 

In  summary,  there  is  no  simple  explanation  for  the  calibration  shifts  observed,  and 
it  is  beyond  the  scope  of  the  program  to  investigate  absolute  calibration  anomalies  of  the 
small  sampling  of  meters  tested.  More  importantly,  we  should  note  that  the  purpose  of  the 
laboratory  tests- is  to  provide  as-installed  calibrations  of  the  sensors  for  use  in  the 
field  tests. 


3.2.2  Field  Test 

The  pipe  assembly  was  buried  7 meters  east  of  the  laboratory  to  a center-line  depth 
of  1.5  meters  and  backfilled  with  clean,  dry  mortar  sand.  As  shown  in  figure  3.10,  the 
pipe  axis  is  horizontal  with  the  sensors  located  adjacent  to  the  vertical  and  horizontal 
axes.  Type  E thermocouples  in  the  sand  back  fill  measure  the  thermal  gradient  of  the 
earth  near  the  pipe.  Thirty  meters  of  pair-shielded  cable  connect  the  thermocouples  and 
sensors  to  the  instrumentation  in  the  laboratory.  The  heat-flow  meter  signals  were  wired 
directly  into  a multi-range,  multi-point  recorder,  and  the  thermocouple  signals  were 
connected  to  a solid-state  ice  reference,  and  then  to  the  multi-point  recorder. 

Backfilling  of  the  pipe  assembly  was  completed  September  17,  1979,  and  data  acquisi- 
tion began  September  17.  To  provide  a base  for  comparison,  we  have  studied  the  zero 
heating  case  first.  Thus,  all  the  field  test  data  given  in  this  report  are  for  zero  power 
input  to  the  electrical  heater.  Heat  input  is  scheduled  for  January  1980. 

Twenty-eight  days  after  burial,  heat  sensor  D2  failed  with  the  resistance  increasing 
from  25  kilohms  to  several  megohms. 

The  variation  of  the  sensor  signal  with  time  is  given  in  figures  3.11,  3.12,  and  3.13 
for  zero  heater  input.  The  meters  are  wired  such  that  outward  radial  heat  flux  produces  a 
positive  emf.  The  associated  earth  temperatures  are  given  in  figure  3.14  and  the  vertical 
temperature  difference  across  the  pipe  (thermocouples  3 and  5)  is  given  in  figure  3.15. 

The  results  clearly  demonstrate  the  effect  of  the  seasonal  (Fall)  cooling  of  the 
earth's  crust.  The  measured  soil  temperature  gradient  of  6 K/m  is  more  than  100  times  the 
geothermal  gradient  [1].  The  soil  heat  flux,  associated  with  this  relatively  high  soil 
temperature  gradient,  is  seen  by  the  pipe  mounted  heat  flux  meters  as  illustrated  in 
figure  3.11  through  3.13.  This  occurs  because  the  thermal  conductivity  of  the  pipe  wall 
is  much  greater  (approximately  100  times)  than  that  of  the  soil.  Thus,  the  soilheat  flux 
enters  the  bottom  of  the  pipe,  is  conducted  along  the  pipe  wall,  and  exits  the  top. 

The  result  of  this  situation  is  that  a time  varying  heat  flux,  which  is  an  order  of 
magnitude  greater  than  that  which  we  wish  to  detect,  is  superimposed  on  the  normal  pipe 
heat  flux.  One  way  around  this  problem  is  to  algebraically  sum  the  response  of  one  or 
more  pairs  of  sensors  mounted  on  opposite  sides  of  the  pipe.  Since  the  earth  heat  flux 
passes  through  the  opposed  sensors  in  opposite  radial  directions,  the  effect  is  largely 
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Figure  3.11 


Field  response  of  meters  A1  and  A2. 
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Figure  3.12 


Field  response  of  meters  B1  and  B2. 
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Figure  3.13  Field  response  of  meters  Cl  and  C2. 
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Figure  3.14  Soil  temperatures  during  heat  flux  meter  field  evaluation. 


DATE 


Figure  3.15  Vertical  temperature  difference  across  pipe  during  heat  flux  meter  field 
evaluation. 
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canceled.  The  preferred  orientation  for  the  sensors  is  on  the  horizontal  axis,  since  the 
earth  temperature  gradients  and  heat  fluxes  are  much  lower  in  this  direction. 

Figure  3.16  shows  the  results  of  such  signal  addition  for  horizontal  sensors  Cl  and 
C2.  The  variation  of  the  sum  of  the  signals  is  about  20  times  less  than  the  variation  in 
the  individual  signals,  and  it  now  appears  that  we  should  be  able  to  detect  changes  in 
heat  flux  of  one  or  two  times  the  nominal  enclosure  value. 


[1]  H.  S.  Carslaw  and  J.  C.  Jaeger,  Conduction  of  Heat  in  Solids,  Oxford  University 
Press,  London,  1959. 
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Figure  3.16  Effect  of  summing  signals  of  horizontally  opposed  sensors. 


53 


NBS-114A  <REV.  9-78) 


U.S.  DEPT.  OF  COMM. 

BIBLIOGRAPHIC  DATA 
SHEET 


1.  PUBLICATION  OR  REPORT  NO. 

NBSIR  80-1637 


2. Oov^t  Accession  No.  3.  fteciptent's  AcC6$sio»No, 


4.  TITLE  AND  SUBTITLE 

Helium  Research  in  Support  of  Superconducting  Power 
Transmission 

Annual  Report 

(October  1978  - September  1979) 


5.  Publication  Date 

December  1980 


ft.  PerfprmiftgOrga^zation  Code 


7.  AUTHOR(S) 

D.  E.  Daney,  Editor 


8.  Performing  Organ.  Report  No. 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

NATIONAL  BUREAU  OF  STANDARDS 
DEPARTMENT  OF  COMMERCE 
WASHINGTON,  DC  20234 


. 4UmtN0. 


11.  Contract/Grant  No. 

433475-AGS 


12.  SPONSORING  ORGANIZATION  NAME  AND  COMPLETE  ADDRESS  rstreer.  City,  state.  ziP) 

Brookhaven  National  Laboratory 
Power  Transmission  Project 
Upton,  NY  11973 


13.  Type  of  Report  & Period  Covered 

Annual , 10/78  - 9/79 


14,  i^nsorSng  Ag«icy  Code 


15.  SUPPLEMENTARY  NOTES 


I I Document  describes  a computer  program;  SF-185,  FIPS  Software  Summary,  is  attached. 


16.  ABSTRACT  (A  200-woTd  or  less  factual  summary  of  most  significant  information.  If  document  includes  a significant  bibliography  or 
literature  survey,  mention  it  here.) 

During  FY  79,  the  NBS  Thermophysical  Properties  Division  program  of  research  for 
superconducting  power  transmission  line  (SPTL)  development  focused  on  three  tasks: 

1)  Analytical,  numerical,  and  experimental  modeling  of  counter-flow  SPTL  cool-down. 

2)  Measurement  of  radial  permeability  of  a SPTL  cable  to  gas  flow. 

3)  Experimental  evaluation  of  thermal  flux  meters  as  a possible  technique  for  deter- 
mining enclosure  heat  leak. 

We  have  completed  the  first  phase  of  our  investigation  of  SPTL  cool-down  in  which  the 
aim  was  to  develop  a fundamental  understanding  of  the  counter-flow  method  of  cool-dowr 
The  excellent  agreement  of  our  analytical,  numerical,  and  experimental  results  gives 
us  confidence  that  we  can  accurately  model  the  cool-down  of  full  scale  SPTL's. 

The  cable  permeability  measurements  indicate  that  any  significant  rupture  of  the  lead 
gas  barrier  of  the  cable  will  lead  to  an  unacceptably  high  leak  rate  from  the  pres- 
surized core. 

During  FY  80  we  shall  concentrate  on  developing  cool-down  strategies  for  full  scale 
SPTL's  using  the  detailed  computer  code  previously  developed  at  NBS.  Field  evaluatior 
of  thermal  flux  meters  will  also  continue. 


17.  KEY  WORDS  (six  to  twelve  entries;  alphabetical  order;  capitalize  only  the  first  tetter  of  the  first  key  word  unless  a proper  name; 
sep^ated  by  semicolons) 


d by  semicolons) 

cool-down,  cable  permeability,  counter-flow  cool-down,  heat  flux  meters, 
liquid  helium,  permeability  of  composites,  superconducting  power  transmission. 


18.  AVAILABILITY  ^ Unlimited 

I I For  Official  Distribution.  Do  Not  Release  to  NTIS 

[3  Order  From  Sup.  of  Doc.,  U.S.  Government  Printing  Office,  Washington,  DC 
20402 

[3  Order  From  National  Technical  Information  Service  (NTIS),  Springfield, 

VA.  22161 


19.  SECURITY  CLASS 
(THIS  REPORT) 

UNCLASSIFIED 


20.  SECURITY  CLASS 
(THIS  PAGE) 


UNCLASSIFIED 


21.  NO.  OF 
PRINTED  PAGES 


58 


22.  Price 
$7.00 


YY  U.S.  GOVERNMENT  PRINTING  OFFICE:  1980-778-159/173  Region  No.  8 


USCOMM-DC 


s. 


